Phase Separation of Binary Systems 



Tian Ma 

Department of Mathematics, Sichuan University, Chengdu, P. R. China 



oo 

o 

O 

(N 

C 

o 



■4— » 

CO 



I 

C 

O 

o 



> 

On 
OO 
(N 

O 

OO 
O 

> 

X 

S-H 



Shouhong Wan^] 

Department of Mathematics, Indiana University, Bloomington, IN 47405^ 
(Dated: June 10, 2008) 

In this Letter, three physical predictions on the phase separation of binary systems are derived 
based on a dynamic transition theory developed recently by the authors. First, the order of phase 
transitions is precisely determined by the sign of a parameter Kd (or a nondimensional parameter 
K) such that if Kd > 0, the transition is first-order with latent heat and if Kd < 0, the transition 
is second-order. Second, a theoretical transition diagram is derived, leading in particular to a 
prediction that there is only second-order transition for molar fraction near 1/2. This is different 
from the prediction made by the classical transition diagram. Third, a critical length scale L c d is 
derived such that no phase separation occurs at any temperature if the length of the container is 
smaller than the critical length scale. 
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Materials compounded by two components A and B, 
such as binary alloys, binary solutions and polymers, are 
called binary systems. Sufficient cooling of a binary sys- 
tem may lead to phase separations, i.e., at the critical 
temperature, the concentrations of both components A 
and B with homogeneous distribution undergo changes, 
leading to heterogeneous spatial distributions. The main 
objective of this Letter is to precisely describe the phase 
separation mechanism and to make a few physical pre- 
dictions. 

Cahn-Hilliard Equation. Let ua and ub be the 

concentrations of components A and B respectively, then 
Ub = 1 — ua- In a homogeneous state, ub = ub is 
a constant. We take u to be the concentration density 
deviation u = ub — Ub- The Cahn-Hilliard free energy is 
given by 



F(u) = F 



/(«) 



dx, 



(1) 



where 



f(u) = a\u 2 + a 2 u 3 + a^u 4 



The same results in this article can be derived in the 
same fashion, and for simplicity, we take this form of / 
as given here. Then the classical Cahn-Hilliard equation 
is as follows: 



du 
di 



= -kA 2 u + A^u 1 + b 2 u 2 + b 3 u 3 }, 



u(x, t)dx = 0, 



(2) 
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supplemented with the Neumann boundary condition: 
du dAu 



dn dn 



o on an, 



(3) 



where n = n^ =1 (0, Lfc) C M 3 is a rectangular domain. 
We note that the more general domain case can be stud- 
ied as well. 

To derive the nondimensional form of equation, let 



A 



l 2 b x 



72 



U = UqU , 

2 b 2 u _ l 2 b 3 u 2 

k ' 73 " k ' 



where I is a given length, uo = ub is the constant con- 
centration of B, and 73 > 0. Then the equation ^ can 
be rewritten as follows (omitting the primes) 



du 

at 



= -A 2 u - XAu + A(7 2 u 2 + 73M 3 ), 
u(x, t)dx = 0, 



u(x, 0) = ip. 



(5) 



Criteria of separation order. Each 3D rectangular 
domain is one of the following two cases: 

Case I: L = L x > L 3 Mj < 2, 3, 

Case II: L = L\ = L 2 > £3 or L\ = L 2 = L3. 

We define a nondimensional parameter: 

2L 2 

2 72 ~ 73 f° r Case I, 
K I (6) 



9tt 2 ' 
26L 2 
27tt 2 



7 2 — 73 for Case II. 



2 



which, by Q 
parameter: 



is equivalent to the following dimensional 



2£|6| 
9tt 2 fc 

27?r 2 fc 



for Case I, 



for Case II. 



(7) 



where Ld = L ■ I is the dimensional length scale. 

By theorems proved in [3], the order of transitions is 
determined by the sign of this parameter K or Kd as fol- 
lows, and we have readily derived the following physical 
predictions: 

Physical Conclusion I: The order of phase separation 
is completely determined by the sign of the nondimen- 
sional parameter Kd as follows: 

(1) If Kd < 0, the separation is second order and the 
dynamic behavior of the Cahn-Hilliard system is as 
shown in Figure^ 

(2) If Kd > 0, the separation is first order transition 
with latent heat. In particular, there are two critical 
temperature T* > T c such that if the temperature 
T > T* , the system is in the homogeneous state, 
when T* > T > T c , the system is in metastable 
state accompanied with hysteresis corresponding to 
saddle-node bifurcation, and when T < T c , the sys- 
tem is under phase separation state. In addition, 
the critical temperatures are functions of uq and L: 
T* = T*(u ,L),T c = T c (u ,L). See Figure^ 
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FIG. 2: For fixed uo and L, the transition for the case where 
Kd > is first order separation with latent heat and with 
hysteresis: Ui and represent separation states, and «o is 



the homogeneous state. In this case, for T c < T < T* , all 
states uo, Ui, u% are metastable states. For T < T c , uo is 
unstable, and Uf and are stable states. 



stable region 




FIG. 3: Typical phase diagram from classical thermodynamic 
theory with shadowed being the metastable region . 



theory (see Reichl [5]), 62 and 63 can be expressed in 
two explicit formulas. Disregarding the term |Vu| 2 , the 
molar Gibbs free energy takes the following form 



FIG. 1: The state uo ~ ub is stable if T c < T, and the state 
Mo is unstable, Uf and U2 are stable if T < T c . 

This is in agreement with part of the classical phase 
diagram from the classical thermodynamic theory given 
in Figure [3jsee, among others, Reichl [8], Novick-Cohen 
and Segal [7J , and Langer [T]. However, as we shall see 
below, our result shows that near uq — 1/2, there is no 
metastable region; see Figure [6| 

Transition diagram. We now examine the order of 
separation in terms of the length scale Ld and mol frac- 
tion uq. For this purpose, according to the Hildebrand 



/ =//^(l — u) + jJbBU + RT(\ — u) ln(l — u) 

+ RTu\nu + au{\ -it), (8) 



where /ia, (J-b are the chemical potential of A and B re- 
spectively, R the molar gas constant, a > the measure 
of repel action between A and B. Therefore, the coeffi- 
cients b 2 and 6 3 are given by 



62 = 



63 = 



D d 3 f(u Q ) 
3! du 3 
D d 4 f{uy) 
4! du A 



2ltn - 1 



6ug(l - 
1 — 3uq 



uo)' 
+ 3u, 



-DRT, 



(9) 



12ug(l-«o) J 



DRT, 



3 



where D is the diffusion coefficient, ft is easy to see that 

f 

2' 
f 

2' 

6 3 > VO < u Q < 1 . 



= if uq = 
^ if u + 



ft is clear that the above formulas for 62 and 63 based 
on the Hildebrand theory fail near uq = 0, 1. However, 
the physically relevant case is away from these two end 
points of mq, and then we have: 



b 2 = 



IQDRT 
3 

ADRT 



o(l). 



(10) 



Then solving Kg = gives a critical (dimensional) length 
scale Ld- 



37TV k V63 

3iry/3k y/b^ 
26 l&al 



for Case f, 



for Case ff. 



(11) 



By (10l and (ff ), we have 



3v 3kn 



L„ 



8^2DRT c \u - || 

9^77 



8-y/26Z?i?T c |u 



0(f) for Case f, 



0(1) for Case II, 



(12) 



where T c is the critical temperature as given in Physical 
Conclusion f. From this formula, we derive the transition 
diagram given by Figure |4j and consequently, we derive 
a theoretical phase diagram given in Figure [5] In partic- 
ular, we have shown the following physical conclusions: 

Physical Conclusion If. 

(1) For a fixed length scale L — L' , there are numbers 
x\ < \ < x 2 such that the transition is second- 
order if the molar fraction X\ < uq < X%, and the 
transition is first-order if Uq > X2 or uq < X\. 

(2) The phase diagram Figure^ is for this fixed length 
scale L' . The points x\ and x-i are the two molar 
concentrations where there is no metastable region 
and no hysteresis phenomena for x\ < uq < x^. In 
other words, 

T*(u ) = T c (u ) for x t < u < x 2 . 

TL-ph&se diagram. We now derive the length and 
temperature phase diagram. For this purpose, we con- 
sider the linear eigenvalue problem for the Cahn-Hilliard 
equation as follows: 



- A 2 u - XAu 
du dAu 
dn dn 



(3u, 
on dn. 



(13) 




FIG. 4: Transition diagram: region II is the first order transi- 
tion region with latent heat, and region I is the second order 
transition region. 



metastable region 




FIG. 5: Derived theoretical diagram for a fixed length scale 
L' with the shadowed region being the metastable region. 



The first eigenvalue is given by 

^2 / „2 



A 



7T / TT 



By (10), we have 

Dd 2 f(u ) 



DRT 



a 
2' 



2 du 2 2u (l — Uq) 

The critical parameter curve equation (3\ — is given by 
tin(l — uq) ( hrr 2 



T r =- 



RD 

Up(l - Uq) 

RD 



2l 2 L 2 

kTT 2 

2^1 



(14) 



Using this formula and the theorems in [4], we derive the 
TL phase diagram given by Figure [6] and the following 
physical conclusions: 

Physical Conclusion III. For a given molar fraction 
< uq < 1, there is a critical (dimensional) length 



T c — 



kir 2 
~2a 



such that the following hold true: 

(1) For < L c d , there is no phase separation for any 
temperature. 



4 




FIG. 6: TL phase diagram. 



(2) For Ld > L c d , phase separation occurs at the critical 
temperature T — T c given by (H). 



Summary. Based on a dynamic transition theory de- 
veloped recently by the authors [21 |3j [6], a system- 
atic mathematical analysis is made for the Cahn-Hilliard 
equation modeling phase separation of binary systems 
3]. Based on this rigorous analysis, we are able to make 
three physical predictions on the phase separation of bi- 
nary systems: 

First, the order of phase transitions is precisely de- 



termined by the sign of a parameter K e i (or a nondimen- 
sional parameter K) such that if Kd > 0, the transition is 
first-order with latent heat and if Kd < 0, the transition 
is second-order. This parameter Kd is explicitly given in 
terms of the system properties and the geometry of the 
container. 

Second, a theoretical transition diagram is derived, 
leading in particular to a prediction that there is only 
second-order transition for molar fraction near 1/2. This 
is different from the prediction made by the classical tran- 
sition diagram. 

Third, a critical length scale L c d is derived such that 
no phase separation even occurs at any temperature if the 
length scale of the container is smaller than the critical 
length scale. The transition temperature T c is precisely 
given as well for the length scale is larger than the critical 
scale. 

Finally, our theory fully reveals the transition dy- 
namics. This is the advantage of using the dynamic 
classification scheme as proposed in [3j |6] , where the 
transitions are classified as Type-I, Type-II and Type- 
Ill. Also, we would like to mention that our results are 
derived for rectangular domains, and more general do- 
main case can be studied using the dynamic transition 
theory as well, and other transition types such as the 
mixed transition may occur; see jl]. 
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